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Abstract. The asymmetric time dependence and various statistical properties of 
polarity reversals of the Earth's magnetic field are utilized to infer some of the 
most essential parameters of the geodynamo, among them the effective (turbulent) 
magnetic diffusivity, the degree of supercriticality, and the relative strength of the 
periodic forcing which is believed to result from the Milankovic cycle of the Earth's 
orbit eccentricity. A time-stepped spherically symmetric a^-dynamo model is used as 
the kernel of an inverse problem solver in form of a downhill simplex method which 
converges to solutions that yield a stunning correspondence with palcomagnctic data. 
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1. Introduction 

The hydromagnetic dynamo in the Earth's outer core converts gravitational and thermal 
energy into magnetic energy |T] . One of the most impressive feature of the geomagnetic 
field is the irregular occurrence of polarity reversals. Averaged over the last few million 
years the mean rate of reversals is approximately 4-5 per Myr, although the last reversal 
occurred approximately 780000 years ago. At least two, but very likely three [3l H] 
superchrons have been identified as "quiet" periods of some tens of millions of years 
showing no reversal at all. 

Knowledge on ancient magnetic field data is mainly obtained from paleomagnetic 
measurements of permanent magnetization from (frozen) lava and sedimentary rocks. 
However, appropriate paleomagnetic sites are rare, unevenly distributed across the 
Earth's surface and, furthermore, actual dating methods allow only a rather rough time- 
resolution. Hence only few reversal characteristics can be evaluated as robust [2]. One 
of the commonly accepted features of reversals is their pronounced temporal asymmetry 
with the initial decay of the dipole being much slower than the subsequent recreation 
of the dipole with opposite polarity [5]. 

Recent numerical simulations have been successful in reproducing not only the 
dominance of the axial dipole and the spectrum of the geomagnetic field, but also the 
irregular occurrence of polarity reversals [SI [71 [8]. Polarity reversals were also observed 
in one [9] of the recent liquid sodium dynamo experiments which have flourished during 
the last decade [TO] . 

It is important to note, however, that neither in simulations nor in experiments it 
is possible to accommodate all dimensionless parameters of the geodynamo [Il],[T2], and 
many of them are not even well known [13j . 

A complementary way to acquire knowledge about the geodynamo is to use available 
magnetic field data for constraining the source of dynamo action. Most famous among 
those attempts is the frozen-fiux approximation [T^ [15] which allows to infer the 
tangential flow at the core-mantle-boundary from secular variation of the geomagnetic 
field. Going beyond this frozen-fiux approximation and trying to infer properties of the 
geodynamo in the depth of the Earth's core has been less successful so far. Such trials to 
"look inside the dynamo" by utilizing spectral data worked nicely for simplified dynamo 
models [161 [IZl [H] but are hardly applicable for real world dynamos. 

With this sobering experience in mind, in the present paper we undertake a 
rather uncommon attempt to constrain some of the most significant parameters of 
the geodynamo by various characteristics of paleomagnetic reversal records. Most 
prominent among those characteristics is the above mentioned temporal asymmetry 
of reversals. Two further characteristics reflect some sort of ordering in the otherwise 
irregular reversal sequence. The first one is the clustering property of reversals which was 
discovered only recently [T9l [20] . Clustering of reversals is manifested in an enhanced 
probability of a consecutive reversal shortly after a first reversal has occurred, which 
results in a deviation from the Poisson distribution that would hold for uncorrelated 
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events. The second one is the appearance of a ~100 kyr periodicity in the distribution 
of the residence times (RTD) between reversals which is beheved to resuh from the 
Milankovic cycle of the Earth's orbit eccentricity [2T1 [22]. 

Based on these three input features we will examine in this paper a simplifying 
a^-model of the geodynamo for which we estimate the degree of supercriticality, the 
noise level, the relative strength of the periodic forcing, and the effective (turbulent) 
magnetic diffusivity of the Earth's outer core. Actually, a few dependencies on individual 
parameters were already published in preceeding papers. In [23] the dependence of 
typical time scales on the supercriticality of the dynamo was studied, in [20] some 
dependencies of the clustering property on the supercriticality and the noise level 
were shown, and in [2l] the influence of the diffusion time scale on the RTD between 
reversals was touched upon. What is new in the present paper is that we take all three 
reversal features together and try to infer from them some essential parameters of the 
geodynamo. 

We will start with a presentation of the forward dynamo problem for which we 
will use a rather simple, spherically symmetric mean-field dynamo of the type. This 
simple model had turned out to be quite helpful for understanding the basic principle 
of the reversal process as a noise-induced relaxation-oscillation in the vicinity of an 
exceptional point of the spectrum of the non-selfadjoint dynamo operator [231 ISSl [26] . 
This exceptional point, at which two real eigenvalues coalesce and continue as a complex 
conjugated pair of eigenvalues, is associated with a nearby local maximum of the growth 
rate situated at a slightly lower magnetic Reynolds number. It is the negative slope of 
the growth rate curve between this local maximum and the exceptional point that makes 
stationary dynamos vulnerable to noise. Then, the instantaneous eigenvalue is driven 
towards the exceptional point and beyond into the oscillatory branch where the sign 
change of the dipole polarity happens. 

After having delineated the simplified mean-field dynamo model we will present the 
solution method for the inverse problem and the main results. 

The paper concludes with a summary and a speculation on the possible 
consequences of our findings for the general understanding and the numerical simulations 
of the geodynamo. 

2. The forward problem 

Before tackling the inverse problem by evaluating the reversal characteristics of many 
different solutions of the time evolution equation for the magnetic field we have to 
delineate in the present section the forward problem, although this has been described 
already in [2^ . 

The governing equation for the mean magnetic field B is the induction equation 
without any mean flows (v = 0) under the influence of a helical turbulence parameter 
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a ^ ^ 

^ = V X (aB) + — V^B . (1) 
OT noa 

This equation results from pre-Maxwell's equations when the source of the magnetic 
field generation is supposed not to be a large scale velocity v but some turbulence which 
comprises helical parts. Since the magnetic field is divergence-free we can decompose it 
into a poloidal and a toroidal parts according to 

B = -V X (r X VS) - r X VT . (2) 

In spherical geometry, the two defining scalars S and T are easily expanded in spherical 
harmonics of degree / and order m. In the following we will assume a to be spherically 
symmetric, being well aware of the fact that this grave simplification does not apply to 
the Earth's outer core. The great advantage of this simplification is that the induction 
equation decouples into pairs of partial differential equations for each degree / and order 
m, 

dsi Id' /(/ + !) ^ , 



dt^ _ld_ 
Ot r dr 



— {rti) ~ a{r,T) — {rsi) 



-[ti- a{r,T)si] . (4) 



^2 



where we have already used dimensionless parameters in which the radius r is measured 
in units of the radius of the Earth's outer core, i?, the time r in units of the diffusion 
time Trf := fiQaR'^, and a in units of (fioaR)"^. The boundary conditions are: 
dsi/dr\r=i + (/ + l)si(l) = ti{l) = 0. In the following we will consider only the dipole 
field with / = 1. 

Due to the presupposed spherical symmetry of a there is no preferred direction of 
the magnetic field axis, hence the order m of the spherical harmonics does not show 
up in the equation system (3-4). This absence of a preferred direction of the dominant 
dipole could be considered a significant weakness of our model. Strictly speaking, our 
restriction to the axial dipole (the mode with / = 1, m = 0) can only be justified if 
some additional symmetry breaking mechanism is supposed to work. This is simply 
the prize we have to pay for the great advantage of remaining in the framework of 
only two coupled partial differential equations (3,4) for si and ti. Certainly, a study of 
more realistic dynamo models would be highly desirable, but with the present computer 
resources it will be difficult to get the good statistical validity that is easily obtainable 
with the large number of reversals resulting from our simple model, let alone the solution 
of an inverse problem as it will be presented in this paper. 

The equation system (3,4), with fixed a(r), would represent a so-called kinematic 
mean-field dynamo model. Below a critical amplitude of a(r), the magnetic eigenfield 
would decay exponentially, above this value it would grow exponentially. In reality, of 
course, the exponential growth of the magnetic field cannot continue indefinitely. After 
having grown to a certain amplitude, the magnetic field attenuates the source of its 
own generation (Lenz's rule). While the precise way of this attenuation is an interesting 
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topic in its own right, we will restrict ourselves to a very simple algebraic "quenching" of 
the kinematic a with the angle averaged magnetic field energy which can be expressed 
in terms of s{r) and t(r). Note again, that this averaging over the spherical angle is 
another simplification that is chosen in order to remain in the framework of a spherically 
symmetric model for which the I and m decoupling of the equation (3,4) remains valid, 
at least formally. 

While the non-linear system of equations that results from (3,4) and the algebraic 
quenching already exhibits a very rich behaviour we will additionally consider the 
infiuence of noise by which a{r) is infiuenced. This noise might be considered as a 
shorthand for fiuctuations of the fiow, changing boundary conditions, and the neglected 
influence of higher magnetic field modes. 

Summarizing the quenching and the noise effect, we model the time dependent 
Q!(r, r) in the form 

a r, r = 



l + E 



+ ii{r) + Ur)r' + Ur)r' + Ur)r\ (5) 

In this equation the noise is considered to have a finite correlation time in which it is 
supposed to be constant. This is equivalent to the following temporal correlation: 

< 6(r)e,(r + ri) > = D\l - \n\/T,) x 6(1 - \r,\/T,)5,, , (6) 

where 6 is the Heaviside step function. In equations (5,6), akin{f) is the kinematic a 
profile, D is the noise intensity, £" is a constant measuring the inverse mean magnetic 
field energy, and Tc is a correlation time of the noise. 

In the following we will motivate the particular choice of the a{r) profile. In former 
papers it was shown that kinematic dynamos of oscillatory character appear only in a 
rather small corridor of a(r) profiles that are characterized by at least one sign change 
along the radius. While this was first shown for the spherically symmetric dynamo 
in [18], a quite similar result was later obtained for a more realistic model in which 
the latitudinal dependence of a is governed by a cos^ dependence [27] . Interestingly, 
such a profiles with one sign change along r were indeed found in simulations of 
magnetoconvection in the Earth's outer core |28j . 

Based on this motivation, we choose for the kinematic a profile in equation (5) the 
particular Taylor expansion 

akin{r) = C ■ (tto + Oiir + a2r'^ + aar^ + a^r'^) (7) 

with «! = as = 0, ^2 = —6 and 04 = 5. The first coefficient, Oq, is chosen close to 1, but 
with two important modifications. First, we add a small parameter 5 which regulates 
the proximity of the kinematic dynamo to oscillatory solutions. Second, ao will also 
incorporate the periodic forcing with the dominant 95 kyr period of the Milankovic cycle, 
being well aware of the fact that the Milankovic cycles contain also other frequencies 
[10]. Hence, taking both effects together we end up with the following ansatz for the 



Inferring basic parameters of the geodynamo from sequences of polarity reversals 



6 



time-dependence of a^: 



a^ir) = 1 + 5 + e cos( 




(8) 



where e parameterizes the strength of the periodic forcing. One could ask why only 
assuming the first Taylor expansion coefficient as time-dependent and not the entire 
expression on the r.h.s. of equation (7). The reason for this is that such a homogeneous 
scaling would have only a weak influence on the dynamo, since it results simply in a 
stronger quenching of akinij) ending up (approximately) at the same quenched a{r) 
profile as before. In contrast to this, a selective change of only one Taylor expansion 
coefficient in equation (7) changes the shape of the a(r) profile, which is much more 
effective for changing the reversal probability. Going back to the real geodynamo, it is 
also not very likely that the Milankovic cycle of the Earth's orbit eccentricity changes 
the dynamo source homogeneously. 

The equation system (3-5), with the concretization (6-8), is time-stepped by means 
of a standard Adams-Bashforth method with radial grid spacing of 0.02 and time step 
length of 2 X 10~^. The correlation time has always been set to 0.005 -Tc which would 
correspond to 1 kyr in case that the diffusion time is set to 200 kyr. The resulting 
time series show reversal sequences quite similar to those of the geodynamo [231 1261 129]. 
For the sake of simplicity, we define a reversal as the sign change of the poloidal field 
component s(l,r) at the outer radius r = 1. Depending on the precise parameters, 
we get typically some 10^ reversals in 1 day CPU time on a normal workstation. As 
an important characteristic of these sequences we will determine the distribution of 
residence times (RTD) Tr between two subsequent reversals. 

3. The inverse problem 

In this section we will try to determine some of the parameters of the presented dynamo 
model in such a way that the resulting reversal sequences fit the characteristics of 
paleomagnetic data as accurate as possible. 

3.1. The parameters to he determined 

Being well aware of some arbitrariness in the particular choice of free parameters in the 
dynamo model, we have decided to choose five parameters which we believe to represent 
some generic features of every dynamo model. 

The first one, with a comparably clear relevance and interpretation, is the magnetic 
Reynolds number C (that is based, however, not on the velocity but on the helical 
turbulence parameter a). The importance of this parameter comes from the fact that, 
as a measure of the supercriticality of the dynamo, it governs the typical time scale of 
the reversal process. This can be understood as follows: While in the rather "quiet" 
regime, when the dynamo is in one polarity, the value of a is quenched approximately 
to the critical one, amidst a reversal (when the magnetic field energy becomes small) 
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a{r) will get close to the unquenched (kinematic) profile. This means that during this 
short time interval the magnetic field dynamics is dominated by the possibly very large 
instantaneous growth rates (and frequencies) of the kinematic a{r) profile. Roughly 
speaking, the higher C the faster the reversal process happens (see, e.g.. Fig. 8 in |23j). 
In addition to this, C governs also the asymmetry of reversals. To understand this 
effect imagine a kinematic a(r) profile which gives an oscillatory dynamo at the critical 
value of C. In this case, "reversals" are nothing than parts of harmonic oscillations. 
With increasing C, i.e. with increasing supercriticality, these oscillations become more 
and more anharmonic (relaxation oscillations) [261 125]. Very often, then, there is 
another transition value of C beyond which the (noise-free) oscillatory dynamo becomes 
stationary again. However, even in this case reversals can be triggered by noise, and the 
asymmetry will still be governed by C. 

The second parameter to be determined is the noise level, represented by D in 
equation (6). Roughly speaking, the larger D, the more frequent the system will leave 
the stable regime and undergo a reversal process. A quantitative interpretation of D is 
non-trivial, since for our case of finite correlation times the relevant quantity is always 
D / -\/T^, so any interpretation of D is only sensible in combination with the correlation 
time of the noise. 

The interpretation of the third parameter 5 in equation (8) is perhaps the most 
intricate one. Formally, 5 describes a constant shift of the kinematic a profile. In the 
purely kinematic regime, the proximity of the dynamo to an oscillatory solution is very 
sensitive on this parameter. In the highly supercritical regime, this sensitivity is reduced 
but there is still some infiuence of 5 on the probability that a reversal occurs. 

The fourth and the fifth parameter refer to the interplay between magnetic diffusion 
and the periodic forcing due to Milankovic cycle of the Earth's orbit eccentricity. This 
problem was already addressed in [21] where we had found that the appearance of several 
clear-cut maxima in the RTD at multiples of 95 kyr is hardly compatible with a diffusion 
time Td := fiQaR"^ = 227 kyr. The latter would follow from recent estimates [30] of the 
molecular conductivity {a = 4.71 x 10^ (f2 m)~^) of the Fe-Ni-Si alloy at the tremendous 
pressure of the Earth's core. Such a large diffusion time would simply "smear out" the 
various maxima in the RTD which are believed to result from a stochastic resonance 
effect. The interesting point is now that this molecular value of the conductivity could 
possibly be reduced by the turbulent motion of the fluid. This so-called P effect has 
been estimated theoretically in various limiting cases [3l]. There have also been claims 
[32] on the measurement of a few percent (3 effect in a turbulent liquid sodium flow with 
magnetic Reynolds numbers Rm up to 8, and some indications of it have been found 
in the Perm dynamo experiment [33] and the Madison dynamo experiment [M], but a 
clear experimental demonstration of it is still elusive. For the geodynamo with its large 
magnetic Reynolds number and high turbulence level, a reduction of a by a factor 2 or 
so is not out of range. This is the reason why we will treat the diffusion time of the 
dynamo model as the fourth free parameter to be determined. 

At the same time we will keep the period of the forcing flxed to 95 kyr (although 
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the actual forcing contains much more frequencies [10]). As indicated aheady in [23] we 
expect that the solution of the inverse problem will give us a diffusion time comparable 
or slightly smaller than 95 kyr, since otherwise the appearance of several maxima in the 
RTD is hardly explainable. 

It is quite natural to consider the strength of this periodic forcing, e in equation 
(8), as the fifth free parameter. 

3.2. The functionals to he minimized and the inversion method 

Since our dynamo model is quite simple, it cannot be expected to reproduce all possible 
reversal features. In particular, spatial distributions of the field during reversals, e.g. 
preferred paths of the virtual geomagnetic pole (VGP), cannot be addressed. However, 
we believe that typical temporal features of reversals, including the shape of individual 
reversals, their clustering properties and the distribution of inter-reversal times, should 
still be reproducible by the model. These three generic characteristics seem to us only 
dependent on the general reversal mechanism and on some basic parameters of the 
geodynamo, five of which were discussed above. 

The first functional to be minimized results from the temporal shape of individual 
reversals which we have adopted from figure 4 of [5] . Actually, the curve " Real" in our 
figure 1 shows the time evolution of the virtual axial dipole moment (VADM) averaged 
over the last five reversals, between -53 kyr and 17 kyr. The first functional Fg^ape 
reflects the deviation of the simulated curves from this "real" curve: 



Here, the subscript real refers to the averaged paleomagnetic measurement, and num 
to the numerically obtained ones. Note that for the latter we have taken an average 
over 100 reversal events, since we thought that an average over only five reversals would 
result in artificial variations with no physical significance. At has been set to 1 kyr, in 
accordance with the resolution in [5]. 

With the second functional we intend to accommodate the clustering property of 
reversal events which had been identified by Carbone et al. in 2005 [19], and analyzed 
in more detail in a follow-up paper by Sorriso-Valvo et al. [20]. The authors had defined 
a quantity which is able to detect clusters or voids even for Poisson processes with time 
dependent rate parameters. Consider the zth event in the sequence of reversals, and 
denote the shorter of the preceding or following inter-reversal intervals (or residence 
times) by 5t, i.e. 



17 



"shape 

(VADM,,„K«At) - VADM„„„(^At))2 



(9) 



i=-53 



5t = min{ti+i -ti]ti- 
Then define the quantity h{5t., At) by 

h{6t, At) = 26t/{26t + At) 



(10) 




wherein 



At — — tj+i if St — — ti 



(12) 
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At = - ti_2 if 5t = ti- ti_i . (13) 

One can easily imagine that in the case of clustering it is rather likely that close 
to the short time interval 5t there will be another short time interval At, hence the 
denominator in equation (10) will acquire a small value and the quantity h{5t, At) will 
be comparably large. Evidently, h takes on values between and 1, and we can either 
consider the probability distribution p{h), or the so-called "surviving function" P{h > 
H) that h is larger than a certain value H. For a Poisson process, with time dependent 
rate parameters, it was shown that the latter probability is Ppoisson{h > H) = 1 — H 
[35] . The actual curve that results from paleomagnetic measurements was shown in 
figure 2 [19], and is also given in figure 2 of the present paper. Based on this, the second 
functional to be considered is 

1000 

Fcluster = E (Pih > H.Ual - P{h > H,)^^j' . (14) 
i=l 

with Hi = i- 0.001. 

While the asymmetry of reversals is rather well understood in terms of the field 
dynamics in the vicinity of an exceptional point of the spectrum, the physical reason of 
the clustering property is less obvious. In [26l [20] we had speculated that the appearance 
of clusters, which can be visualized as a "devil's staircase" of reversal events [36], is a 
typical sign of "punctuated equilibrium" [37] and self-organized criticality [38]. In this 
respect it is also interesting to note the tendency of the clustering property to increase 
with the degree of super-criticality |20j . 

The third functional refers to the distribution function of inter-reversal times (or 
residence times). The appearance of certain periods in the reversal sequence had been 
controversially discussed for a long time [391 SOI SI] • A breakthrough was achieved by 
Consolini and DeMichelis who analyzed not the Fourier spectrum but the probability 
distribution of inter- reversal times [21]. In this distribution they observed a clear 
sequence of several maxima at multiples of approximately 95 kyr. The physical reason 
for these maxima was identified as a stochastic resonance with the Milankovic cycle of 
the Earth's orbit eccentricity [22]. The data from [21] are shown again in figure 3, and 
the corresponding functional is 

90 

4 = 1 

with Ti = i ■ 1 kyr. 

Obviously, having defined the three functionals, there still remains an ambiguity in 
the choice of the relative weights for them. We will test some reasonable choices of the 
a-priori errors a shapes dusters and artd in the total functional 

Ftotal = ^ shapeFshape + clusterFcluster + ^rtdFrtd (16) 

and check afterwards the correspondence of the resulting curves with the paleomagnetic 
ones. 

The very inversion is carried out by using a standard downhill simplex method 
taken from " Numerical Recipes" [12] . Since we have 5 parameters to be determined, we 
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Table 1. Parameters of the three versions resuhing from the downhill simplex 
inversion. 





C 


5 


Ccrit 


C / Ccrit 


D 




£ 


#/Myr 


Version 1 


57.0 


-0.03 


13.0 


4.4 


6.3 


86.8 kyr 


0.133 


2.05 


Version 2 


96.8 


-0.16 


8.99 


10.8 


8.0 


63.7 kyr 


0.123 


3.07 


Version 3 


147.8 


-0.18 


8.73 


16.9 


7.8 


55.1 kyr 


0.094 


2.91 



use a simplex with 6 points. For each of those points in parameter space we solve the 
forward problem for 20000 diffusion times which typically gives a few thousand reversals. 
From these reversals the functional (16) is computed, and based on this evaluation the 
usual steps of the downhill simplex method are performed, until an appropriate stopping 
criterion is reached. While an individual run takes a few hours on a normal PC, the 
solution of the inverse problem takes about one week. 

3.3. Results 

In the following we will present three solutions which result from choosing different 
weights of the three individuals functionals in (16). 

Table 1 shows the obtained parameters of the three versions. The second raw gives 
the parameter C, the third raw the shift parameter 6. In the fourth raw the critical 
value of the dynamo for the a(r)-profile with ao = I + S, ai = = 0, a2 = —6 and 
a4 = 5 (see Eqs. 5-7) is shown, which differs from version to version due to the variation 
of 6. The fifth raw gives the degree of supercriticality, i.e. C/ Ccrit, where Ccrit prescribes 
the amplitude of the critical a at which the onset of dynamo action occurs. The raws 6, 
7 and 8 show the noise level D, the (effective) diffusion time T^, and the strength of the 
periodic forcing e, while the last raw gives the resulting number of reversals per Myr for 
each version. 

Figures 1, 2 and 3 show the curves for the temporal dependence of the VADM, 
for the probability P{h > H), and for the RTD between reversals, compared with the 
corresponding curves from the paleomagnetic data. 

Obviously, version 2 gives the best correspondence with all three reversal 
characteristics. Compared to this "optimal" version 2, in version 1 Crtd was increased 
by a factor 2 which results in a deteriorated correspondence with the real RTD curve in 
figure 3. In version 3, aghape was increased by a factor 2 (compared to version 2), which 
leads to a deteriorated correspondence with the real VADM curve in figure 1. Evidently, 
the large value of C and the small Td provide reversals that are faster than the really 
observed ones. 

4. Conclusions 

With version 2, we have obtained the best reproduction of the paleomagnetic input 
data for a 10 times supercritical dynamo, a relative strength of the periodic forcing 
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Figure 1. Temporal evolution of the virtual axial dipole moment (VADM) during 
reversals. Comparison of the average over the last five real reversal data (taken from 
[5]) with the average over 100 reversals for the three versions resulting from different 
solutions of the inverse problem. 




Figure 2. Surviving function P{h > H) for the real data (adopted from [HlllS]) and 
for the three versions. The deviation from the straight line which would appear for a 
Poisson process indicates a significant clustering. 



of some 15 per cent, and an effective magnetic diffusion time of approximately 65 kyr 
which is by a factor 3.5 smaller than the value that would result from the molecular 
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T [Myr] 

Figure 3. Residence time distribution (RTD) for the real data (taken from [2T|) and 
for tlie three versions. 

conductivity. The latter is perhaps the most important result of the inversion, for 
the following reason: The conductivity enhancing effect of turbulence is a subject of 
ongoing debate. Only very few estimations exists for the turbulent diffusivities within 
the Earth's outer core. However, their values are of extraordinary importance for the key 
parameters of the geodynamo (e.g. Ekman number or the magnetic Prandtl number) 
as well as for the estimation of turbulent transport properties and the destructive 
influence of turbulence on magnetic field generation through turbulent field diffusion. 
The molecular values of the dimensionless parameters are extremely small and cannot 
be realized in simulations so that usually enhanced values are applied that shall resemble 
the effective (i.e. turbulent) quantities. A more detailed information on the turbulent 
diffusivities therefore delivers important information on the significance of the parameter 
space accessible in global MHD simulations for the geodynamo. Furthermore, reality is 
more complicated, as it is very likely that a conductivity reduction due to turbulence 
would be anisotropic, because the small scale turbulence in a fast rotating object like the 
Earth is subject to preferred directions parallel to the rotation axis (and also along the 
dominating field component) so that the convection cells (that determine the turbulent 
transport of physical quantities) are oriented parallel to the rotation axis and elongated 
along the magnetic field. It is well known that an anisotropic conductivity could have a 
tremendous effect on the selection between equatorial and axial dipole solution H5] . 
Roughly speaking, axially aligned rotating columns tend to decrease the conductivity 
for horizontal currents less than the conductivity for vertical currents, and this effect 
leads to a preference of the axial dipole compared to the equatorial dipole. This could be 
an extremely important point since the conclusions of many geodynamo related papers 
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rely on an isotropic conductivity when looking for criteria for the selection of axial or 
equatorial dipoles. 

The stunning agreement of paleomagnetic and numerical reversal characteristics as 
shown in figures 1 , 2 and 3 gives support to our hypothesis that reversals are indeed noise 
triggered relaxation oscillations in the vicinity of an exceptional point of the spectrum 
of the dynamo operator. In this respect it is important to note that in particular the 
time asymmetry and the clustering property are intrinsic and robust properties of the 
model that appear for very wide regions of parameter (if not for all). By solving the 
inverse problem we have not "produced" them, but have only fine-tuned the dynamo 
parameters to fit optimally the paleomagnetic data. 

We have carefully tried not to over-interpret our simple model by focusing only on 
those parameters to be determined, and those functionals to be minimized, that refer 
to the temporal properties of reversal sequences, and not to any spatial features. This 
makes us optimistic that the results will prove robust when inversions of this kind will 
later be repeated using more realistic dynamo models. Given that one downhill simplex 
run for our simple model takes already one week, one can imagine that corresponding 
runs with better dynamo models will lead to significant computational costs. 
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